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ABSTRACT 

Recent measurements of the Sunyaev-Zel'dovich (SZ) angular power spectrum from the South Pole 
Telescope (SPT) and the Atacama Cosmology Telescope (ACT) demonstrate the uTiportance of un- 
derstanding baryon physics when using the SZ power spectrum to constrain cosmology. This is 
challenging since roughly half of the SZ power at £=3000 is from low-mass systems with lO^'^/i"^ 
Mq < M500 < 1-5 X \Q^^h~^ Mq, which are more difficult to study than systems of higher mass. 
We present a study of the thermal pressure content for a sample of local galaxy groups from Sun et 
al. (2009). The group Ysph,500 — -^500 relation agrees with the one for clusters derived by Arnaud et 
al. (2010). The group median pressure profile also agrees with the universal pressure profile for clusters 
derived by Arnaud et al. (2010). With this in mind, we briefly discuss several ways to alleviate the 
tension between the measured low SZ power and the predictions from SZ templates. 



Subject headings: cosmology: observations — galaxies: clusters: general 
— cosmic background radiation 



X-rays: galaxies: clusters 



1. INTRODUCTION 

Galaxy groups are not scaled-down versions of rich clus- 
ters following simple self-similar relations (e.g., Ponman 
et al. 2003; Voit 2005). They are systems where the role 
of complex baryon physics (e.g., cooling, galactic winds, 
and AGN feedback) begins to dominate over gravity. The 
effects of these baryon processes are not large but still 
significant in massive clusters, and therefore need to be 
calibrated if we want to further improve the cosmological 
constraints from clusters. Since the role of these processes 
is less pronounced in massive clusters, it is easier to study 
and understand them by observing groups. The impor- 
tance of galaxy groups for cosmology has also been well 
demonstrated by the recent measurements of the Sunyaev- 
Zel'dovich (SZ) angular power spectrum. SPT measured 
a value for the SZ power at £ = 3000 (scales of ~ 4') that 
was lower than most prior predictions by at least a factor 
of two (Lueker et al. 2010). A measurement of the SZ 
power spectrum by ACT is consistent with these results 
(Das et al. 2010; Dunkley et al. 2010). The SZ angular 
power spectrum is a sensitive probe both of cosmological 
parameters and of the hot gas content of galaxy clusters 
and groups (e.g. Komatsu & Seljak 2002). Regarding the 
latter, it opens a new window into low-mass, high-redshift 
systems as about half of the SZ power at £=3000 comes 
from halos with IQ^^h''^ Mq < M^qq < 1.5 x lO^^h'^ Mq 
and z > 0.5 (e.g., Trac et al. 2010). While the exami- 
nation of the thermal pressure content in z > 0.5 groups 
is a challenge to current X-ray telescopes with typical ex- 
posures, such work can be done for local groups. Group 
pressure profiles have received little attention to date and 
the existing samples are small (e.g., Mahdavi et al. 2005; 
Finoguenov et al. 2007). In this letter, we present the 
pressure profiles of hot gas in 43 local galaxy groups from 
the Sun et al. (2009, S09 hereafter) sample. We assume 



17m=0.24, Oa=0.76 and Hq = 73 km s^^ Mpc'^ 

2. THE GROUP SAMPLE AND DATA ANALYSIS 

The group sample and the Chandra data analysis have 
been discussed in S09. There are 43 groups at z=0.01 - 
0.12 (a median z of 0.033), all with intracluster medium 
(ICM) properties derived to at least r25oo- Twenty-three 
have masses measured to ~ ''soo- The mass range is 
M500 = lO^^/j-i M© - 10"/i"^ Mq. As an archival sam- 
ple, there is not a well defined selection function. The 
sample does include some X-ray faint groups, usually with 
strong radio AGNs at the center. 

The analysis in S09 was based on Chandra CALDB 
3.4.3. Since then, there have been major Chandra calibra- 
tion releases on the on-axis effective area and the ACIS 
low-energy contamination models. To check the impact of 
calibration changes on the S09 results, we examined spec- 
tra from 40 regions from 18 observations which we studied 
in S09. The temperatures of these regions range from 0.7 
to 3.0 keV and the abundances range from 0.15 to 1.5 so- 
lar. The observation dates are from 2000 to 2006, and both 
ACIS-I and ACIS-S observations are examined. CIAO 4.3 
with CALDB 4.4.1 and XSPEC 12.6 were used. The tem- 
perature decrease with the new calibration is less than 1% 
on average. This is not surprising (e.g., Nevalainen et al. 
2010) as temperatures at this range are mainly determined 
by the centroid of the iron-L hump. The decrease of the 
ACIS effective area below 5 keV causes the normalization 
of the spectral model to increase, on average, by 9.9%, 
which is independent of the temperature. This implies an 
average 4.9% increase in gas density. We do not include 
this small change in this work. 

While small ACIS calibration uncertainties may still re- 
main after the release of CALDB 4.4.1, the results for 
groups should not be affected much. The Chandra flux 
in the 0.5 - 2.0 keV band agrees with the XMM-Newton 
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flux to better than 4% (e.g., Nevalaincn ct al. 2010). The 
existing temperature bias is smaU for low-temperature gas 
as long as the iron-L hump is significant. Even if temper- 
atures of low-temperature gas are still biased high by e.g., 
10%, the l^ph,500 ~ ^500 relation (Fig. 1) should not be 
affected much as M500 and ysph,500 will both be lower, by 
~ 15% and ~ 20% respectively, producing a small change 
in the relation that is insignificant compared to the large 
statistical uncertainties on M500 (a median la uncertainty 
of 18%). 

3. THE INTEGRATED PRESSURE CONTENT 

We first compare the scaling relation between mass and 
the volume-integrated Compton parameter for the S09 
groups with the Arnaud et al. (2010, AlO hereafter) re- 
sults. The AlO sample is from XMM- Newton observations 
of 31 Representative \MM- Newton Cluster Structure Sur- 
vey (REXCESS) clusters with M500 = 7 x IQ^^hr^ M© 
- 6 X 10^^h~^ Mq. The spherically integrated quantity, 
i^sph,500 (defined in Equ. 14 in AlO), is derived within rsoo- 
The S09 Fsph.500 ^ ^500 relation agrees well with AlO's 
(Fig. 1). This can be expected from the good agreement 
of the M500 — ^,500 relation between the two works. The 
difference between Ix. ioo and Vsph.soo is the ratio between 
the spectroscopic temperature and the gas-mass-weighted 
temperature within rsoo {Tx vs. Tmg,5oo)- The Tx in S09 
is measured between 0.15 rsoo and rsoo (we call it Tsoo), 
while the Tx in AlO is measured between 0.15 rsoo and 
0.75 rsoo. For the S09 sample, Tmg.soo / ^500 = 0.98±0.04 
and Tsoo = (0.95±0.02) rx(0.15 rsoo - 0.75 rsoo)- Arnaud 
et al. (2007) quoted 0.94-0.97 for the latter ratio. Thus, 
rmg,50o / Tx(0.15 rsoo - 0.75 rsoo) ~ 0.93 for the S09 sam- 
ple, which agrees with the AlO result of 0.924±0.004. Be- 
cause of this consistency and the agreement between the 
M500 - lx,500 relations, we expect the Fsph.soo - Afsoo re- 
lations from S09 and AlO to agree. We also examined the 
twenty tier 3 and 4 groups in S09 (where gas properties 
are only derived to 45% - 72% of rsoo)- The Kiph,2500 - 
Tsoo relation for the tier 3 and 4 groups agrees with the 
relation for the tier 1 and 2 groups. Overall, we conclude 
that for the S09 groups, the derived l^ph,500 is 1.05±0.25 
times Fsph.500 predicted from the AlO relation. We also 
plot the predicted ygph.soo ^ Afsoo relations from recent 
SZ templates (Schgal et al. 2010a; Battaglia et al. 2010; 
Shaw et al. 2010; Trac et al. 2010) in Fig. 1. The Seh- 
gal et al. (2010a) template used the ICM model by Bode 
et al. (2009), which was calibrated with the Vikhlinin et 
al. (2006) and S09 gas fraction relations. Vikhlinin ct al. 
(2009) included six more clusters with on average lower gas 
fractions within rsoo than the eight clusters in Vikhlinin et 
al. (2006) . As the Msoo — ^,soo relation from Vikhlinin et 
al. (2009) agrees well with that in AlO, it is not surprising 
that the Sehgal et al. (2010a) template has ~ 9% higher 
normalization than the AlO result for clusters. 

4. THE RADIAL PRESSURE PROFILES 

We also examine the radial pressure profiles of the S09 
groups. AlO derived a universal pressure profile of the 
ICM by removing the mass dependence. If the ICM scal- 
ing relations are self-similar, the mass dependence is M^/^ 
(Equ. 5 of AlO, Psoo, also see Nagai et al. 2007). Since 
deviations from self-similarity exist, AlO defined a term 



(Equ. 7 and 8 of AlO, which we call ^adjust) to further 
remove the mass dependence in addition to PsoO) where 
^adjust « M"oo"^ and a{x) = 0.22/(l-f2;3), x = 2r/r5oo. 
While the form of Padjust ean be examined with the S09 
groups, the large uncertainties in Afsoo for the S09 groups 
do not allow good constraints on this adjustment factor. 
We derived the P/ Psoo and the P/Psoo /Padjust profiles for 
all 43 groups. Uncertainties are estimated from those in 
the temperature and density profiles. In Fig. 2, we show 
the median and the la scatter of both profiles. To ac- 
count for uncertainties in pressure profiles, 1000 Monte 
Carlo simulations were run to examine the uncertainties 
on the median profiles, which is small as shown in Fig. 2. 
The la scatter for P/Psoo/-Padjust is 26% - 40% at 0.2 rsoo 
- 0.8 rsoo (compared to less than 30% scatter for the AlO 
clusters) and is consistent with a log-normal form. Thus, 
the radial pressure profiles of the S09 groups agree well 
with the imiversal pressure profile defined by AlO. This is 
consistent with the good agreement between the two works 
on the mass proxies and entropy scalings (S09; Pratt et al. 
2010). 

5. DISCUSSIONS AND CONCLUSIONS 

The results above suggest that the thermal pressure of 
local galaxy groups from the S09 sample is consistent with 
the extrapolation from the AlO results, although statis- 
tical errors and scatter are still large. Interestingly, re- 
cent measurements of the SZ angular power spectrum are 
at least a factor of two lower than prior expectations at 
£ = 3000 (Lueker et al. 2010). The thermal SZ power 
spectrum scales roughly as the square of the thermal SZ 
flux. Given the results presented above, we briefly discuss 
several possibilities that may alleviate this tension. 

X-ray selection bias: The S09 sample is an X-ray 
archival sample and the REXCESS sample is an X-ray- 
luminosity-selected sample. Both samples can be different 
from mass-selected samples. The Chandra archive may be 
biased to systems with bright cores, while X-ray under- 
luminous groups and clusters may exist (e.g., Rasmussen 
et al. 2006; Popesso et al. 2007). However, the F^ph - M 
relation at rsoo and beyond is less affected by the presence 
of X-ray bright regions (e.g., a large cool core) than the 
Lx- M relation. For the S09 sample, 12% - 68% of the 
X-ray flux (a median of ~ 34%) is from within 0.15rsoO) 
while such regions only contribute ^ 5% to K;ph.50o for 
Msoo = 10"/?^! - lO^^h'^ Mq halos, assuming the AlO 
pressure profile. The contribution may be even smaller 
than 5% for X-ray under-luminous systems as their gas 
cores are fainter than those of the REXCESS clusters used 
to derive the AlO profile. One main conclusion of the S09 
work is that the gas content of groups is comparable to 
that of clusters at r > r2S00j at least for the S09 sample. 
If we combine the — Tsoo relations from Vikhlinin et al. 
(2009), S09 and REXCESS, the trend of slope flattening 
with increasing radius is significant, with an almost con- 
stant density at rsoo from groups to clusters. This trend is 
consistent with the scenario that much of the low-entropy 
gas in low-mass systems has been ejected to large radii by 
strong feedback (e.g., McCarthy et al. 2010). However, it 
remains to see whether this result applies to mass-selected 
samples. One way to test this is to examine scaling rela- 
tions from non- X-ray-selected samples. This kind of work 
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has been done on optically selected samples by stacking 
the ROSAT all-sky survey data (Dai et al. 2010; Rykoff 
et al. 2008). Besides the systematic uncertainties with 
stacking, the ROSAT temperatures from stacking are of- 
ten biased (Rykoff ct al. 2008) and contamination to those 
samples, especially at the low-mass end, can be severe. 

Pressure contribution at r > rsoo: The total SZ flux 
is more sensitive to the gas in cluster outskirts (r > rsoo) 
than the total X-ray flux. Few direct X-ray constraints ex- 
ist at such large radii, especially for groups. Although the 
contribution from r > r5oo to l^ph is significant (e.g., 40% 
- 70% increase by integrating to 2r5oo for the AlO profile), 
the contribution to the SZ power spectrum at ^—3000 as- 
suming the AlO profile is smaller, only about 20% from 
r > rsoo regions. So an overestimate of the thermal pres- 
sure from the AlO model only at r > rsoo could overpredict 
the SZ power spectrum by at most 20% at ^=3000. 

Dynamical state of the ICM: The SZ signal mea- 
sures the total thermal energy of electrons. However, the 
potential energy of halos may not be fully converted into 
thermal energy of electrons because of e.g., recent merg- 
ers and weak viscosity of the ICM (e.g., Lau et al. 2009; 
Burns ct al. 2010). Galaxies can also contribute to the 
non-thermal pressure support by e.g., injected magnetic 
fields and cosmic rays. The non- thermal pressure support 
may also cause the ICM to be clumpy. For a clumpy ICM, 
the SZ signal predicted from the X-ray data will be biased 
high. All thcisc cffcxrts may have a dependence on mass (or 
the ICM temperature), and the evolution of these effects 
with redshift may not be self-similar. 

The impact of the non-thermal pressure support on the 
SZ power spectrum has been discussed by Shaw ct al. 
(2010), who examined models with a radial dependence 
of the non-thermal pressure. Trac et al. (2010) examined 
a model with 20% non-thermal pressure for all clusters 
and groups at all masses and rcdshifts. However, both 
models do not predict the moan value measured for the 
SZ power spectrum at ^ = 3000 by Luekcr ct al. (2010), 
being high by about la. Interestingly, SZ observations 
suggest the latter model predicts too little SZ flux for very 
massive clusters (Schgal ct al. 2010b). As for dumpiness, 
the good agreement between the measured SZ radial pro- 
file and the prediction from X-ray data for individual clus- 
ters (e.g., Plagge et al 2010; Komatsu et al. 2010; Sehgal 
et al. 2010b) suggests that dumpiness should be weak 
for massive clusters. However, one can imagine a mass 
dependence for dumpiness, as both heat conduction and 
dynamic viscosity can be much weaker in groups than in 
clusters. 

If the hydrostatic equilibrium mass under-estimates the 
true mass by 20% from M500 = 10" h'^ Mq - 10" h"i 
Mq (e.g., Nagai et al. 2007), 54ph,500 will increase by 4% - 
7%, if the universal pressure profile from AlO is assumed. 
Therefore, the normalization of the l^ph.soo ^ Afsoo rela- 
tion will decrease by ^ 24%, which roughly translates to 
~ 42% decrease on the predicted SZ power spectrum. If 
the mass bias is larger for low-mass halos, the relation will 
be steeper and the decrease will be larger. However, it 
is a big challenge to constrain the effects of non-thermal 
pressure support in groups and clusters, especially its de- 
pendence with mass and redshift. For groups, robust mass 
measurements that do not assume hydrostatic equilibrium 



are required. Two promising methods are stacking of the 
lensing data (e.g., Leauthaud et al. 2010) and caustics 
(Rines & Diaferio 2010). Future X-ray microcalorimeter 
observations (e.g., by Astro-H) may also constrain the ICM 
turbulence directly. 

Evolution of the ICM properties: While local 
groups are discussed in this paper, most of the SZ power 
at £=3000 from groups is from systems at z > 0.5. Evolu- 
tion of the ICM properties is poorly constrained for poor 
clusters and groups. Recent results on z > 0.5 groups sug- 
gest that the evolution of the Lx scalings is not weaker 
than the self-similar prediction (e.g., Jeltema et al. 2009; 
Leauthaud et al. 2010), but the statistical uncertainties 
arc large. Proper understanding of the evolution of the 
Lx scalings requires a good understanding of the selection 
function (e.g. Pacaud et al. 2007). Better constraints on 
the evolution of the low-mass end of the Lx scaling re- 
lations should be achieved with more XMM-A^ewfon and 
Chandra data on larger samples with well-defined selec- 
tion functions. Of course, the evolution of the Lx scaling 
relations is not equal to the evolution of the y — M re- 
lation. More factors, e.g., the evolution of the cool core 
fraction and the gas distribution, need to be accounted 
for. Alternatively, deep SZ observations of high-z groups, 
either individually or by stacking, can directly constrain 
the evolution of the Y — M relation. 

Contamination from radio and infrared galaxies: 
Both radio and infrared galaxies could potentially fill in 
SZ decrements at 150 GHz. While the contamination from 
radio galaxies should be small (see discussion in Sehgal 
et al 2010a), the contamination from dusty star-forming 
(infrared) galaxies is less clear. In Lucker et al. (2010), 
the signal from infrared galaxies was removed from maps 
at 150 GHz by subtracting maps at 220 GHz after fitting 
for a weighting factor. If all infrared sources have the same 
spectral index of a = 3.6, then this should effectively re- 
move infrared contamination from the 150 GHz maps. If 
some infrared sources have a shallower slope (e.g., a = 2.6 
as in Knox et al. 2004), then residual contamination will 
remain. However, a more recent analysis by Shirokoff et 
al. (2010) suggests that even a large correlation between 
infrared galaxies and groups/clusters would not increase 
the 95% CL upper limit on the thermal SZ power spec- 
trum to the level that it is consistent with predictions 
prior to Lueker et al (2010). 

This work shows that the local groups from the S09 
sample follow the extrapolation of the pressure scaling re- 
lations from AlO. More data are required to reduce the 
statistical errors and more importantly explore the system- 
atic uncertainties discussed above. Regarding the low SZ 
power measured by recent experiments, we suggest some 
astrophysical possibilities that may alleviate the apparent 
tension between models and measurements. Understand- 
ing the SZ power spectrum will provide important insights 
into both baryon physics and cosmology. 

Wc thank S. Allen, N. Battagha, P. Bode, G. Holder, 
J. Hughes, D. Nagai, J. Ostriker, L. Shaw, J. Sievers and 
H. Trac for helpful discussions. We thank N. Battaglia, 
L. Shaw and H. Trac for providing their best-fits of the 
Y — M relations. M.S. and C.S. were supported in 
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Fig. 1. — Upper panel: the Kjph,500 ~ -'^SOO relation for the S09 groups (black points) and the AlO clusters (red points). The best-fit 
relation for the S09 groups (the black solid line) is derived from the BCES orthogonal regression method with bootstrap resampling (Akritas & 

Bershady 1996), and is given by £;(^)"^/^Kjph,500 = 10*^ (A^SOO / 3x10^^ Mq)° h~^''^ Mpc^, where a = 1.75±0.09 and /3 = -4.77±0.09. 
The black dotted lines show the Icr error (22%). The red dotted line shows the AlO best-fit relation. The hydrostatic equilibrium (HSE) 
mass values for the REXCESS clusters are not published so the Afsoo — ix.soo relation (Equ. 2 in AlO) was used to derive M500, which is 
the reason for the small errors and scatter of red points as yx,500 ^sph,500 s-re well correlated. Good agreement between the S09 and 
AlO results can be seen. We also plot the best-fit relations from recent SZ templates; the blue dashed line is from Sehgal et al. (2010a) 
(a = 1.80, fS = 4.72), the green dashed line is for the AGN feedback simulations at z = by Battaglia et al. (2010) (a = 1.73, /3 = 4.81), 
the magenta dashed line is from Shaw et al. (2010) at z = 0.05 (a = 1.81, /3 = 4.76 for the HSE mass, a = 1.80, /3 = 4.85 for the true mass) 
and the cyan dashed line is for the nonthermal20 model by Trac et al. (2010) (a = 1.83, /3 = 4.86). We emphasize that M500 from the X-ray 
data is the HSE mass which may be smaller than the true M500. The models by Shaw et al. (2010) and Trac et al. (2010) assume about 
20% non-thermal pressure support. If plotted with the HSE mass, these two lines will shift ~ 12% to the left for Afsoo (or ~ 23% higher for 
^sph,50o)- Lower panel: the ratios between the SZ templates and the AlO best-fit (the same color code as in the upper panel), while the black 
solid line shows the ratio between the S09 and the AlO best-fits. 
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Fig. 2. — Normalized pressure profiles of 43 groups from S09. P500 is the normalization factor to remove the mass dependenee if the 
ICM scaling relations are self-similar (Equ. 5 of AlO), while Padjust is the factor defined in AlO (Equ. 7 and 8) to further remove the mass 
dependence as the observed ICM relations deviate from self-similarity. The blue solid and dotted lines show the median and la scatter of the 
group pressure profiles. The green dotted lines show the Icr uncertainties of the median. The three red lines in the left panel are the universal 
pressure profile from AlO for M500 = 10^^ h^^ Mq, 10^'* h^^ Mq and 10^^ h^^ Mq, from the top to the bottom respectively. The median 
-'^500 of the S09 groups is ~ 7 X 10^^ h^^ Mq. The red line in the right panel is the universal pressure profile from AlO after removing the 
mass dependence. The median pressure profile for the S09 groups agrees well with the AlO profile. The dashed line shows the position of 
J'2500- 



